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Abstract 


Efforts  are  made  to  relate  the  depolarized  light  scattering  spectral  power  density  to 
the  dynamic  shear  compliance  spectrum  of  a  viscoelastic  liquid.  Generalized  Langevin 
equations  are  used  to  describe  the  dynamics  and  coupling  of  reorientation  and  translational 
motions.  It  is  shown  that  in  the  presence  of  the  coupling,  the  depolarized  scattering 
spectrum  consists  of  two  components,  one  is  associated  with  the  reorientational  motion,  and 
the  other  with  the  viscoelastic  relaxation  of  the  solution.  Methods  to  obtain  the  dynamic 
shear  compliance  from  the  depolarized  light  scattering  spectrum  are  described.  It  is  also 
shown  that  in  the  limit  of  very  slow  reorientational  motion,  the  depolarized  light  scattering 
spectrum  is  proportional  to  the  imaginary  part  of  the  dynamic  compliance  spectrum,  a  result 
predicted  by  the  Pockels  photoelastic  law. 
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L  Introduction 


Relating  the  dynamic  light  scattering  spectral  density  to  the  mechanical  relaxation 
spectrum  is  an  important  area  of  research.  Making  the  correlation  brings  about  the 
unification  of  information  that  is  currently  obtained  separately  by  the  two  techniques,  hence 
enabling  one  to  seek  the  desired  physical  data  with  one  technique  when  the  other  becomes 
inconvenient  or  failed.  The  dynamic  light  scattering  spectral  density  in  a  bulk  liquid  occurs 
as  a  result  of  the  thermally  induced  modulation  of  the  dielectric  constant  tensor  due  to  local 
density  and  orientational  fluctuations  of  molecules  inside  the  scattering  volume.  If  the  bulk 
liquid  is  a  polymer,  the  local  density  and  orientational  fluctuations  are  associated  mainly 
with  the  motion  of  chain  segments.  On  the  other  hand,  mechanical  relaxation  occurs  as  a 
result  of  segmental  motion  and  as  such  it  is  expected  to  be  related  to  dynamic  light 
scattering.  Both  the  local  density  fluctuations  and  the  orientational  fluctuation  contribute 
to  the  polarized  ( W)  component  of  the  scattered  light,  whereas  the  depolarized  (VH)  light 
scattering  intensity  is  associated  with  the  orientational  fluctuation. 

In  an  one  component  system,  it  has  been  shown  that  the  dynamics  of  density 
fluctuation  is  equivalent  to  the  relaxation  of  longitudinal  compliance.^  Consequently, 
dynamic  light  scattering  serves  as  an  alternative  technique  to  measure  the  longitudinal 
compliance  spectrum,  which  has  generally  been  obtained  by  a  mechanical  relaxation 
spectrometer. 

In  the  binary  polymer  solution  consisting  of  a  polymer  and  a  small  molecular  weight 
solvent,  owing  to  mutual  diffusion  arising  from  osmotic  pressure  fluctuations,  the  isotropic 
component  of  the  light  scattering  spectrum  is  dominated  by  the  concentration  fluctuation. 
However,  because  of  the  mixing  of  the  concentration  fluctuation  and  viscoelasticity,  the 
viscoelasticity  that  affects  the  mechanical  properties  of  the  polymer  solution  is  also  expected 
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to  contribute  to  the  dynamic  light  scattering  spectrum.  The  effect  of  viscoelasticity  on  the 
dynamic  light  scattering  spectrum  of  the  binary  polymer  solution  was  recently  analyzed  by 
this  author.^^ 

In  depolarized  light  scattering  (DLS),  the  scattered  intensity  arises  from  the 
anisotropy  fluctuations  of  the  dielectric  constant  tensor  as  a  result  of  molecular 
reorientation.'^  Because  of  the  technological  importance,  attempts  have  been  made  to  relate 
the  DLS  spectrum  to  the  d>uamic  mechanical  properties.  Hess  attempted  to  relate  the  DLS 
spectrum  to  the  mechanical  relaxation  spectrum  by  using  the  Pockels  photo-elastic  law,  in 
which  the  off  diagonal  part  of  the  optical  impermeability  tensor  is  set  proportional  to  the 
shear  strain  tensor.^  Since  the  Pockels  photo-elastic  law  completely  neglected  the 
contribution  due  to  molecular  reorientation,  it  is  only  applicable  to  the  description  of  the 
deformation  of  the  condensed  medium  owing  to  acoustic  waves  at  high  frequency  at  which 
molecular  reorientation  fails  to  follow.  In  viscous  molecular  liquids,  it  is  well  known  that 
reorientation  of  molecules  makes  a  dominant  contribution  to  the  DLS  intensity.^’^  In  the 
viscoelastic  region,  both  the  molecular  reorientation  and  overdamped  shear  waves  contribute 
to  the  DLS  spectral  density  and  they  can  be  probed  by  the  dynamic  light  scattering 
technique.^’^  For  viscoelastic  liquids,  such  as  polymer  melts  and  polymer  solutions, 
measurements  of  the  dynamic  shear  modulus  and  compliance  are  of  technological 
importance.^®  It  is  thus  of  interest  to  find  a  method  with  which  information  concerning  the 
dynamic  shear  compliance  and  shear  modulus  data  can  be  extracted  by  the  optical  means 
using  the  DLS  spectrum.  Focusing  on  this  objective,  we  present  in  this  paper  the  result  of 
a  theoretical  effort  which  succeeds  in  correlating  the  DLS  spectrum  with  the  dynamic  shear 
compliance  of  a  viscoelastic  liquid.  The  experimental  verification  of  the  theoretical  result 
wll  be  the  subject  of  our  future  publication.  In  section  II  we  provide  the  theoretical 
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background  needed  for  carrying  out  the  analysis.  In  section  III  we  consider  the  coupled 
motions  of  the  off-diagonal  part  of  the  polarizability  density  and  translation  in  the 
siscoelastic  liquid.  In  section  FV  we  show  how  the  dynamic  shear  compliance  can  be 
obtained  from  the  depolarized  light  scattering  spectrum. 


II.  Theoretical  Background 

Consider  a  beam  of  electromagnetic  radiation  with  wave  vector  kj  and  frequency  oj 
propagating  through  a  medium  characterized  by  an  anisotropic  dielectric  constant  tensor^ 
consisting  of  an  average  part  e  and  a  fluctuating  part  6^ . 

The  spectral  power  density  of  the  scattered  light  I^j  (g,<«>)  is  given  by"^ 

e-‘“‘  <a^(q,t)  a*(q)> 


where  jg  is  the  scattering  vector  given  by£  =  k^  -  kj ;  k^  is  the  amplitude  of  k^,  the  wave 
vector  of  the  scattering  light,  and  is  equal  to  27rn/XQ.  Here  n  is  the  refractive  index  of  the 
medium  which  is  related  to  €  by  n  =  Aq  is  the  wavelength  of  the  incident  light  in 

vacuum,  r  is  the  radial  distance  from  the  scattering  center  to  the  detector.  The  quantity 
S  Cjj  (g,o))  is  the  projection  of  the  fluctuating  part  of  the  dielectric  constant  tensor  onto  the 
incident  light  polarization  indicated  by  and  the  scattered  light  polarization  indicated  by 

A 

°s- 

o)  is  related  to  the  polarizabiUty  density  (g,  w)  by 


5e^(g,o)  =  4n 


(2) 


If  the  molecules  are  weakly  coupled  electronically,  then  the  polarizability  density  tenser 

Page  5 


component  a^j  (g,t)  given  above  can  be  written  as: 


=  E  “2^(0exp{iq*r^(t)} 

where  is  the  local  field  conected  polarizability  tensor  component  of  molecule  k 
located  at  positions  Both  and  depend  on  time  because  of  molecular 
reorientation  and  translational  diffusion.  OK'ing  to  Eq.  (3),  measurements  of  the  spectral 
power  density  I^j  (q.o)  involving  with  different  polarization  vectors  of  the  incident  and 
scattered  light  beams  will  yield  information  about  molecular  motion  in  a  dense  fluid. 

Using  the  scattering  geometry  with  the  incident  light  propagating  along  the  x-axis 
and  polarized  along  the  y-axis  of  the  laboratory  coordinate  system,  it  can  be  shown,  for  a 
system  consisting  of  identical  optically  anisotropic  molecules  dissolved  in  an  optically 
isotropic  solvent,  that  the  linear  polarizability  density  involved  in  depolarized  (VH)  light 
scattering  at  wave  vector  jq  has  two  components:^ 

«VH  =  “yz  (9/2)  +  (q,t)  sin  (0/2)  (4) 

where  0  is  the  scattering  angle,  and  q  is  the  amplitude  of  the  scattering  vector  given  by 

47tnsm(6/^  expressions  for  and  for  molecules  of  arbitrary  shape  are  most 


conveniently  given  in  terms  of  rotational  matrices.^’^^  However,  for  optically  symmetric  top 
molecules,  (g,t)  and  (g,t)  are  given  by: 

^  Sny(j)n^(j)  e'a  'Ij 

j 
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(6) 


“xy(9-0  =  ^2i^(j)Qy(j)e‘^  '-j 

j 

where  ^  =  aj  -  aj_,  aj  and  a_L  being  the  components  of  the  local  field  conected  molecular 
polarizability  parallel  and  perpendicular  to  the  SNtninetry  axis,  respectively.  The  position 

A 

coordinate  and  orientation  of  the  symmetry  axis  of  molecule  j  are  labeled  byjj  and  n  (j), 
respectively.  The  quantity  n^(j)  is  the  projection  of  the  unit  vector  representing  the 

A 

symmetry  axis  of  molecule  j  onto  the  laboratory  a-axis,  a  =  x,  y,  z.  Both  and  n  (j)  are 
time  dependent,  but  we  do  not  explicitly  write  out  the  time  dependence  of  jj  and  n^  (j)  in 
Eqs  (5)  and  (6)  for  brevity. 

Using  Eq.  (4)  one  can  readily  show  that  the  depolarized  (VH)  spectrum  for  a 
rotationally  isotropic  medium,  is  given  by: 

m 

IvH  (q>“)  =  (V2  tr)  J  dt  e*‘“^  {  (g,t)  (9)>  cos2(e/2) 

+  <“xy  (3’^)  “xy’(3)>  sin  2(0/2)}  (7) 

where  the  cross  terms,  <a^  (q,t)  a^*(q)>  and  <a^  (g,t)  do  not  contribute  to 

the  spectral  intensity  because  of  different  reflection  symmetry  for  and  with  respect 
to  the  yz  plane.  It  should  be  noted  that  Eq.  (9)  is  valid  for  molecules  with  arbitrary 
geometry,  not  being  restricted  to  symmetry  top  molecules. 

III.  Coupling  of  Rotation  and  Translational  Motions 
In  order  to  calculate  the  time  conelation  functions  associated  with  the  VH  spectrum, 
we  need  to  consider  the  types  of  molecular  motion  that  affect  the  time  dependence  of 
(q,t)  and  (q,t).  Andersen  and  Pecora^  considered  this  problem  in  1971  in  the  theory 
of  the  depolarized  scattering  spectrum  from  overdamped  shear  wave  in  a  viscous  molecular 
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liquid.^^  They  used  the  symmetry  consideration  and  concluded  that  the  conserved 
hydrodynamic  variables,  such  as  number  density  and  energy  density  will  not  couple  to  both 
and  However,  because  of  identical  reflection  symmetry  associated  with  the  variables 
to  Py  and  the  y  component  of  the  linear  momentum  density  Py  will  couple  to  Oy^.  From 
this,  they  then  constructed  a  microscopic  theory  to  describe  the  dynamics  of  the  coupled 
variables  using  the  technique  of  generalized  hydrod>'namics.^^ 

To  relate  the  viscoelastic  effect  to  the  depolarized  light  scattering  spectrum  of  the 
fluid,  we  generalize  in  this  paper  the  two-variable  theory  of  Andersen-Pecora  by  including 
the  memory  effect  in  the  relaxation  kernel.  We  use  a  vector  variable  A(g,t)  consisting  of 
a  set  of  two  dynamic  variables,  one  of  which  is  Oy^  (q,t)  and  the  other  is  Py(q,t).  As  given  - 
by  Eq.  (2),  the  crj^(q,t)  variable  is  related  to  the  cooperative  orientational  motion  of 
molecules  within  the  scattering  volume: 


A(q,t) 


iP/q.t)J 


(8) 


where  Py(q,t)  is  the  y-component  of  the  linear  momentum  density  given  by: 

Py(q,t)  =  S  ry(i)eWj  (9) 

j 

Here  rty^^  is  the  y-component  of  the  linear  momentum  of  molecule  j.  The  time  dependence 
in  TTy^^  and  are  not  explicitly  written  out.  For  an  isotropic  fluid  system,  the  physical 
property  does  not  depend  on  the  direction  of  jg;  only  the  amplitude  is  considered. 

We  need  next  to  obtain  an  equation  for  the  correlation  function  <Oy^(g,t)ay2*(q)> 
in  order  to  calculate  the  depolarized  light  scattering  spectrum.  Having  chosen  the  set  of 
dynamic  variable  A(g,t),  we  can  proceed  to  obtain  the  equations  of  motion  for  and  Py, 
Owing  to  the  fact  that  is  coupled  to  Py  we  shall  obtain  two  coupled  equations.  However, 
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in  contrast  to  the  treatment  of  Andersen-Pecora,  we  keep  the  memory  effect  in  the 
generalized  Langevin  equations.  The  memory  effect  plays  an  important  role  in  determining 
the  viscoelasticity  of  the  liquid.  The  various  essential  steps  needed  for  obtaining  the 
generalized  Langevin  equations  for  and  Py  are  given  in  the  monograph  by  Wang,^^  and 
they  are  briefly  described  in  Appendix  I. 

From  Appendix  I,  the  equations  of  motion  for  and  Py  in  the  small  q  limit  are 
given  by: 


ayx(q.t)  =  -(  dx  Kn(x)ay^(q,t-T)-iq  jf  dTKij(x)Py(q,t-T) 


(10) 


P/q,t)  =iq/^‘<JrK2,(T)ay,(q,t-T)  V /jK^(t)P,(q.t-T) 

where  the  dots  above  and  Py  indicate  taking  the  time  derivative.  The  explicit 
expressions  for  K|j(r)  are  given  in  Appendbc  I.  In  Eqs.  (10)  and  (11),  we  have  neglected  the 
random  force  terms  as  these  terms  vanish  upon  taking  the  ensemble  average  with  respect 
to  the  initial  variables  when  calculating  the  time  correlation  functions.  If  the  various  Kjj(T) 
terms  relax  significantly  faster  than  the  dynamic  variables  Oy^  and  Py,  then  the  Markov 
approximation  can  be  used  to  simplify  Eqs  (10)  and  (11).  If  the  Markov  approximation  is 
assumed,  then  Eqs  (10)  and  (11)  reduce  to: 


ay,(q.t)  =  -ra^(q,t)-iqpPy(q,t) 


(12) 


Py(q,0 


Pm 


(13) 


where  r  is  the  orientational  relaxation  rate  equal  to  the  time  integral  of  Kjj(t);  is  the 


mass  density  and  is  the  shear  viscosity  of  the  liquid  which  is  proportional  to  the  time 
integral  of  K22(t);  m  and  are  rotation-translation  coupling  constants  associated  with  the 
time  integrals  of  Kj2(^')  and  K2i(r),  respectively.  Equations  (12)  and  (13)  are  equivalent 
to  the  equations  used  by  Anderson  and  Pecora  for  the  calculation  of  the  depolarized  light 
scattering  spectrum  of  the  viscous  molecular  liquid.^  Since  we  do  not  assume  the  Markov 
approximation  in  the  present  work,  the  viscoelastic  response  of  the  liquid  is  retained  in  Eqs 
(10)  and  (11). 

To  proceed,  we  need  to  solve  the  coupled  equation  given  by  Eqs.  (10)  and  (11). 
These  are  coupled  differential-integral  equations  of  the  convolution  type  and  can  be  solved 
by  using  the  usual  Laplace  transform  technique. 

Taking  the  Laplace  transform  of  Eqs  (10)  and  (11)  and  after  rearranging,  we  obtain; 


s+q^ 

A  A 

where  cy  (q»s)  and  Py  (q,s)  are  the  Laplace  transform  of  (q,t)  and  Py  (q,t),  respectively. 
They  are  given  by: 


(15a) 

Py(q^)  =  /’Py(q>0e'^dt 

(15b) 

To  obtain  the  correlation  functions,  we  first  invert  Eq.  (14)  to  obtain  the  solution. 
We  next  multiply  ci^(q,s)  by  its  initial  value  a*yz(q),  and  Py(q,s)  by  its  initial  value  Py*(q) 
and  then  take  the  ensemble  averages  for  both  of  them.  The  results  are: 
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(16a) 


Cp(q^)  =  <Py(q^)P;(q)> 


— !— (s +K„)NmkT 
A(q)  " 


(16b) 


where  N  is  the  total  number  of  molecules  and  m  is  the  mass  of  the  molecules.  The  other 
quantities  are  given  by: 

i  (q)  =  (s  t  K„)  (s  ■*  q%)  -q'  Kj^^,  (Ha) 

A  A 

We  can  combine  (s)  and  K21  (s)  and  introduce  a  quantity  C  (s): 

-Kij(s)Kjj(s)  =  {(s)r(s)/p„  (17d) 

The  quantity  C(s)  has  the  same  dimension  of  viscosity  and  will  be  referred  to  as  the  dynamic 
coupling  viscosity  as  its  presence  arises  from  rotation-translational  coupling.  The  djmamic 
coupling  viscosity  C  is  easily  shown  to  be  given  by: 

C(s)  = - ^ -  1  f "  dte"^  <io_(i)  a„(0)>o  p  (18e) 

Vr(s)(kTzJkT  'Jo  yxw  y,v  y  Q  1  v  / 

The  quantity  Xa  appearing  in  Eqs  (17b)  and  (17e)  is  the  static  polarizability 


susceptibility  defined  in  Appendix  I. 


rv.  Depolarized  Light  Scattering  Spectrum  and  Dynamic  Modulus 
We  now  relate  r^s)  to  the  dynamic  shear  (stress)  modulus  G*(<i)).  Using  Eqs  (16a) 
and  (17c)  and  setting  s  =  iw,  we  can  rewrite  the  momentum  density  correlation  function 
given  in  Eq.  (16b)  as: 


Cp(q4w) 


_ NmkT _ 

[io+(q^/pJf|]  +  [q^Cr/pJ/[ico  +  T] 


(18) 


At  zero  frequency,  the  dynamic  coupling  viscosity  C  vanishes,  and  one  obtains  the 
limiting  hydrodynamic  result  for  all  molecular  fluids: 


Urn  Cp(q4(o)  =  (NmkT)/(q^/pJri,  (19) 

o  ->0 


where  is  the  macroscopic  shear  viscosity  equal  to:^ 

q  =  lim  [G"(u))/{o] 
o-O 


(20) 


At  finite  frequency,  the  dynamic  coupling  viscosity  is  finite,  and  from  Eqs.(18)  and  (19)  we 
expect  that 


[q(io))  +  ((io))]  =  G(o)) 


(21) 


A  ^  A  $  _ 

where  G(o)  is  related  to  the  dynamic  shear  modulus  G  by;  io  G  (o)  =  G  (o).  The 
dynamic  modulus  can  be  separated  into  a  real  G'((i))  and  an  imaginary  part  G"(w): 


G*(o)  =  G'(o))  +  iG"(<.>)  *  io)  r“dte'‘“‘G(t)  (22) 

Jo 
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In  the  limit  of  small  frequency,  as  given  by  Eq  (20)  G''(w)/t^  becomes  the 
macroscopic  viscosity  but  G'((i>)/(i)  vanishes.^ 

On  the  account  of  Eq.  (21),  we  can  rewrite  Eq.  (17a)  as: 


and 


A(ia)) 


(io  +r) 


icoC  ' 
io)  +r^ 


(23) 


Cr  (io)  = 


[io)+(qVpJ(G-C)lkTx, 


(io  +r)(iw  +(q^/Pa)G-(qVpJ[iwC/(iw  +r)]) 


(24) 


We  expect  in  general  the  dinamic  coupling  viscosity  C  to  be  much  less  than  the  dynamic 
shear  viscosity  G.  Thus,  we  proceed  by  expanding  Eq.  (24)  in  power  series  of  C/G,  and 
retain  only  the  linear  term  in  C  to  arrive  at: 


cr(io)  = 


1 


cr 


-  (-^) - 


kTx, 


(25) 


For  the  fluid  in  the  viscoelastic  regime,  the  macroscopic  viscosity  is  high  so  that  the 
molecular  reorientation  rate  is  slow.  In  addition,  we  shall  also  neglect  the  dispersion  of  the 
collective  reorientation  rate  r  and  the  dynamic  coupling  viscosity  C;  i.e.  we  shall  assume  that 
r  and  C  are  real  and  independent  of  frequency  in  the  frequency  window  of  the  dynamic  light 
scattering  experiment.  This  assumption  is  justified  experimentally.^ 


We  now  combine  Eq.  (25)  with  the  xy  component  of  the  scattering  spectrum  (c.f.  Eq. 
(7))  and  obtain  for  o  >  r,  the  spectral  power  density  of  the  depolarized  Rayleigh  scattered 


(26) 


=  .  (3l,(X)  eos^l  rJ - ifiC - 

(0^-r^  Pa  2  |-0)2+q2G*(o)/p^ 

where  G*  (o)  is  the  d>T3amic  shear  modulus  given  in  Eq.  (22);  the  symbol  Re  indicates 
taking  the  real  part  of  the  expression  within  the  curvy  brackets. 

One  notes  in  Eq.  (26),  and  more  generally  Eq.  (18),  the  fact  that  the  dynamic  shear 
modulus  enters  the  depolarized  scattering  spectrum  is  due  entirely  to  the  presence  of 
rotation-translational  coupling.  If  the  dynamic  coupling  viscosity  C  vanishes,  then  the 
depolarized  spectrum  becomes  a  simple  Lorentzian,  with  the  half-width  at  the  half-height 
equal  to  r,  the  collective  reorientational  rate.  This  is  consistent  with  the  previous  result.^’^ 
For  the  viscoelastic  liquid,  the  molecular  reorientation  rate  is  slow  and  the  Lorentz  term  is 
very  narrow.  In  addition,  one  also  notes  that  the  second  term  in  Eq.(26)  containing  the 
dynamic  shear  modulus  term  is  proportional  to  cos^  0/2.  Thus,  the  dynamic  shear  modulus 
contribution  is  pronounced  at  small  scattering  angles,  but  it  vanishes  when  the  scattering 
angle  0  equals  180®.  However,  except  for  low  frequencies,  one  notes  that  the  second  term 
is  proportional  to  q  ,  it  also  becomes  small  at  small  scattering  angle.  Hence,  when  the 
factors  of  q^  and  cos^  0/2  are  combined,  one  arrives  at  the  conclusion  that  0  =  90®  is  the 
optimum  scattering  angle  to  observe  the  shear  modulus  contribution  to  the  depolarized 
scattering  spectrum.  However,  if  frequencies  are  low,  the  q^  factor  is  canceled  out  by  the 
q  term  m  the  denominator  and  in  this  situation  the  DLS  spectrum  becomes  independent 
of  the  scattering  angle.  But,  due  to  the  cos^  0/2  factor,  the  intensity  of  the  scattered  light 
associated  with  the  shear  modulus  is  most  pronounced  at  small  scattering  angle. 

Returning  to  Eq.  (26),  we  separate  G*  into  real  and  imaginary  parts.  By  taking  the 
real  part  of  Eq.  (26)  and  after  carrying  out  some  algebra,  we  obtain  the  DLS  spectrum  as: 
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(27) 


f.  ^  r  iBwlv  (q/PmK5 

SvH(q>^)='  ,  (~)Cr(cOS^-)(— ) - ; — — - — ,  -  - 

Pxn  2  o  (-co^+q^G7pj'*(q'G"/pj2 

If  the  depolarized  light  scattering  experiment  is  carried  out  by  using  an  electrom'c 
correlator,  whose  dynamic  range  is  usually  less  than  1  MHz,  then  for  polymer  solutions  with 
a  storage  modulus  greater  than  lO^Pa,  the  conditions  G'/Pm  >  ^  is  valid  over  the 
entire  dynamic  range  of  the  correlator.  In  such  a  circumstance,  the  term  w  m  the 
denominator  of  the  second  term  of  Equation  (27)  can  be  neglected,  and  as  a  result,  Eq.  (27) 
reduces  to: 


Svh(‘i.“)  =  -^^+cr(cos=|)J"M/u 

where  J"(w)  is  the  imaginary  part  of  the  dynamic  shear  compliance,  given  by 


I»(c)  - - 9M -  (29) 

[a'(<j)p  *  [G"(«)]" 

Using  the  definition  of  J'(<i>),^®  we  can  transform  Eq.  (28)  back  to  the  time  domain  as: 

g™(52)  =  «■"  +  (r(cos^|)  ^ -/(()]  (30) 

where  is  the  steady  state  shear  comph’ance,  and  J(t)  is  the  shear  creep  compliance. 
gvH(q.t)  is  the  normalized  VH  polarzability  correlation  function.  Modem  dynamic  light 
scattering  experiments  often  employ  a  digital  electronic  correlator  and  directly  measures  the 
time  correlation  function  gypj  (q,t).  Equation  (30)  is  more  convenient  than  Eq.  (28)  which 
deals  with  the  spectral  power  density. 

Equation  (28)  or  Eq.  (30)  is  an  important  result;  it  shows  that  the  dynamic  shear 
compliance  can  be  obtained  from  the  DLS  spectrum.  Owing  to  the  presence  of  rotation- 
translation  coupling,  the  DLS  spectrum  of  an  optically  anisotropic  viscoelastic  liquid  will 
consist  of  two  components.  One  is  associated  wth  reorientational  motion,  and  the  other 
with  the  viscoelasticity  of  the  liquid.  The  viscoelastic  component  is  proponional  to  the 


imaginary  part  of  the  dynamic  shear  compliance.  This  result  is  quite  analogous  to  that 
previously  obtained  for  the  polarized  scattering  spectrum,  in  which  it  is  shown  that  the 
density  fluctuation  spectrum  as  probed  in  isotropic  scattering  is  proportional  to  the 
imaginary  part  of  the  dynamic  compression  compliance.^  However,  due  to  the  presence  of 
reorientation  motion  which  in  some  cases  contributes  significantly  to  the  DLS  intensity,  the 
depolarized  light  scattering  spectrum  is  proportional  to  the  imagmary  part  of  the  dynamic 
shear  compliance  only  if  the  collective  reorientation  rate  r  is  so  small  that  it  can  be 
neglected.  If  the  reorientational  rate  is  neglected,  then  one  obtains  a  result  identical  to  that 
obtained  by  using  the  photo-elastic  law,^  and  the  factor  Cr  becomes  the  Pockels’ 
photoelastic  constant. 

If  the  reorientational  dynamics  contributes  to  depolarized  light  scattering,  then  the 
reorientational  effect  has  to  be  separated  out  in  order  to  use  depolarized  light  scattering  as 
a  techm'que  to  obtain  the  dynamic  shear  compliance  data.  To  separate  the  reorientational 
contribution  to  the  spectrum,  one  can  first  cany  out  the  DLS  experiment  at  large  scattering 
angle  near  180°.  As  one  notes  in  Eq.  (28),  the  spectrum  at  large  scattering  angle  is  mainly 
due  to  reorientatioiL  After  obtaining  the  angular  independent,  the  reorientational  spectrum 
(or  the  conelation  function,  Eq.  (30))  can  be  subtracted  off  from  the  total  DLS  spectrum 
to  obtain  the  viscoelastic  spectrum. 

An  alternative  method  can  also  be  used.  This  method  allows  one  to  obtain  both  the 
real  and  imaginary  parts  of  the  dynamic  shear  compliance  spectrum.  We  now  describe  this 
alternative  method  of  separating  the  reorientational  contribution  to  obtain  the  compliance 
spectrum. 

A 

From  Eq.  (16a)  we  take  the  inverse  of  (q,s)  and  obtain  by  setting  s  =  io. 
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l(l) 


(31) 


- h -  =io+r+(5^) 

IJ*  dte-^-’Cf  (q,t)  Pm  -<*>^  +  (qVpJG‘ 

To  obtain  Eq.  (31),  we  also  neglect  the  higher  order  terms  in  C  due  to  the  fact  that  C  is  in 

A 

general  much  less  than  G 

We  next  write  the  one  sided  Fourier  transform  of  the  polarizability  time  correlation 
function  appearing  in  Eq.  (31)  as: 

^ /•  die-'”'  cr(q,t)  ”  I(q.<e) -iH(q,e>)  (32) 

Here  I  (q.u)  and  (q.“)  tlie  cosine  and  sine  transforms  of  the  normalized  correlation 
function  (q,t),  respectively.  The  cosine  transform  of  (q,t)  is  proportional  to  the 
yz  component  of  the  depolarized  light  spectrum  given  in  Eq.  (7). 

Using  Eq.  (32)  and  expressing  G*  in  terms  of  real  and  imaginary  parts,  we  obtain 
from  Eq.  (31) 


1 

I(q,0))-iH(q,o) 


ico  +r 


_ io _ 

-o2  +  (-3i)(G'+iG") 

Pm 


(33) 


Equating  the  real  and  imaginary  parts  on  both  sides  of  Eq.  (36)  and  also  neglecting  the 
term,  we  obtain  for  the  real  part 


A(q)  s 


I 


and  for  the  imaginary  part 


r+wcrj"(w) 


(34) 


B(w)  H  =  G>+o>rCJ'w)  (35) 

where  the  frequency  dependence  for  I  and  H  in  Eqs.  (33)  and  (34)  are  not  written  out 
explicitly;  J'  and  J"  are  real  and  imaginary  parts  dynamic  shear  compliance  J*  (w),  given 
by 
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J’(o)  =  JV)-iJ"(w) 

E>ynamic  shear  modulus  G*  (o)  and  dynamic  shear  compliance  J*  (o)  are  related  by  the 
equation^® 

J*(oj)G*(a))  =  1  (37) 

One  notes  that  B((i>)  vanishes  at  o  =  0.  Equations  (34)  and  (35)  can  be  used  to 
separate  the  reorientational  rate  r  from  the  DLS  spectrum.  The  method  consists  of 
computing  the  sine  and  cosine  transforms  of  the  polarizability  correlation  function  that  one 
measures  by  using  a  digital  electronic  correlator.  From  these  one  can  calculate  A(ti))  and 
B((i)).  From  the  plots  of  A  (o)  as  a  frequency  of  q  one  obtains  r,  as  an  intercept.  After 
obtaining  r,  a  one  can  substitute  it  into  Eqs.  (34)  and  (35)  to  calculate  CJ’  and  CJ'  as  a 
function  of  frequency. 

However,  as  one  notes  in  Eq.  (7),  the  depolarized  Rayleigh  spectrum  has  two 
components.  The  component  proportional  to  the  sin^  0/2  term  is  only  associated  with  the 
rotational  dynamics.  The  contribution  of  this  term  must  be  either  subtracted  off  or 
minimized  in  order  to  use  the  results  given  in  Eqs.  (34)  and  (35).  This  can  be  accomplished 
by  carrying  out  the  DLS  experiment  at  a  small  scattering  angle. 

In  conclusion,  we  show  that  it  is  possible  to  use  the  method  of  the  DLS  experiment 
to  obtain  J'  and  J'  or  J(t)  of  a  viscoelastic  liquid.  We  have  used  the  generalized  Langevin 
equations  to  analyze  the  connection  benveen  the  DLS  spectrum  and  dynamic  shear 
compliance.  We  have  shown  that,  owing  to  rotation-translational  coupling,  both 
reorientational  motion  and  viscoelasticity  contribute  to  the  DLS  spectrum.  In  the  limit  of 
negligible  molecular  reorientation,  the  DLS  spectrum  is  shown  to  be  proportional  to  the 
imaginary  part  of  the  dynamic  compliance  spectrum,  a  result  predicted  by  the  photoelastic 
law.  When  both  the  reorientational  motion  and  viscoelasticity  contribute  to  the  DLS 
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intensity,  we  have  provided  methods  for  separating  the  reorientational  component  from  the 
DLS  spectrum  to  obtain  the  dynamic  shear  compliance  data. 
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Appendix  I 


Given  the  coupled  dynamic  variables  (q,t)  and  Py  (q,t)  defined  in  Eqs.  (3)  and 
(9),  we  express  then  as  a  column  matrix  A(q,t). 


A(q.t)  = 


[Pv(q.t). 


(A-1) 


where  in  (q,t)  is  the  yzth  component  of  the  polarizability  density  tensor.  Its  time 
dependence  is  due  to  translation  and  it  describes  the  reorientational  motion  of  a  collection 
of  molecules  inside  the  scattering  volume.  Py  (q,t)  is  the  y  component  of  the  linear 
momentum  density  for  the  collection  of  molecules  inside  the  scattering  volume.  We  follow 
the  techmque  of  generalized  hydrodynamic  theory  and  obtain  the  generalized  Langevin 
equation  of  motion  for  A  (q,t)  as 


-|A(q,t)  =  iQ-A(q,t)-  /‘  K(q.T)-A(q.t-T)dT  +  F(q,t) 
ot  -  JO  " 


(A-2) 


In  Eq.  (A-2),  ^  is  the  frequency  matrix,  and  is  a  null  matrix  for  the  chosen  set  of  dynamic 
variables  given  in  Eq.  (A-l)."^  I^(q.^)  is  the  relaxation  matrix  given  by: 


K(q.T)  =  <F(q.T)F(q)><AA*>- 


(A-3) 


where  F(q,r)  is  the  random  force  which  can  be  calculated  according  to: 
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F(q.t)  =  e‘(‘-P)*(i-p) 


iSfpay,(q) 

F.(q.t) 

iSfP/q)  , 

"  [Fp(q.t) 

1 

A 

0 

0 

1/NmkT 

0  1/NmkT 

where  j  is  the  Liouville  operation  and  p  is  the  projection  operator.  Since 


(A-4) 


(A-5) 


F.(q)  =  i2a„(q)  =  E  [i®*iqK®/nip(r^e‘« 

j 

and  from  the  momentum  conservation  equation,  we  have 

F,(q)  =  ii£P,(q)  =  |p,U  = 

In  the  small  q  limit,  the  stress  tensor  a„  becomes  q-independent,  and  is  given  by 


j 

and  the  initial  values  of  (q)  and  Fp  (q)  become 


(A-6) 


(A-7) 


(A-8) 


F.(<i)=E 


(A-9) 


Fp(q)  =iqo, 


(A-IO) 


Using  Eqs.  (A-4),  (A-5),  (A-9)  and  (A-10),  we  obtain  from  Eq.  (A-3)  the  relaxation 
matrix  j^(q,  r)  in  the  limit  of  small  q  as 


Kji(x)  -iqK,j(T) 
5W'iqK„W  q=K,,(t) 


(A-11) 


where  in  the  limit  of  small  q. 
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Kn(^)=<“y,W«;*(0)Q>/kTX. 

(A-12) 

with 

=  <IE 

j 

(A- 13) 

Ku(t)=<ay,(t)(iOy,)*>/NmkT 

(A-14) 

K2iW=<i‘'yzW«;z(0>/KTx. 

and 

(A-15) 

K^W  =  <o„(l)o;>/NmW  (A-16) 

Here  the  subscript  Q  in  Eq  (A'12)  designates  that  the  time  dependence  is  taken  with 

respect  to  projected  Uouville  space. 

If  we  write  down  from  Eq.  (A-2)  the  equation  for  each  individual  dynamic  variable, 

we  obtain 

■|a^(q,t)  =  -/J<lTKjj(T)ay^(q,t-T)-iq|^‘dTK,2(x)Py(q,t-c) 

■^F.(q,t) 

and 

(A-17) 

|p,(q,t)=iq/JdTK„(t)«^(q.t-T)-,*/jK,,(T)P,(q.t-T) 

which  are  equivalent  to  Eqs.  (13)  and  (14)  provided  that  the  random  force  terms  (q,t) 

and  Fp  (q,t)  are  neglected. 
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